Space-Dependent Step Features: Transient Breakdown of Slow-roll, Homogeneity and Isotropy 

During Inflation 
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A step feature in the inflaton potential can model a transient breakdown of slow-roll inflation. Here we gen- 
eralize the step feature to include space-dependence, allowing it also to model a breakdown of homogeneity and 
isotropy. The space-dependent inflaton potential generates a classical curvature perturbation mode characterized 
by the wavenumber of the step inhomogeneity. For inhomogeneities small compared with the horizon at the step, 
space-dependence has a small effect on the curvature perturbation. Therefore the smoothly oscillating quantum 
power spectrum predicted by the homogeneous step is robust with respect to sub-horizon space-dependence. For 
inhomogeneities equal to or greater than the horizon at the step, the space-dependent classical mode can domi- 
nate, producing a curvature perturbation in which modes of wavenumber determined by the step inhomogeneity 
are superimposed on the oscillating power spectrum. Generation of a space-dependent step feature may there- 
fore provide a mechanism to introduce primordial anisotropy into the curvature perturbation. Space-dependence 
also modifies the quantum fluctuations, in particular via resonance-like features coming from mode coupling to 
amplified superhorizon modes. However these effects are small relative to the classical modes. 

PACS numbers: 98.80.Cq, 04.62.+V 



I. INTRODUCTION 

Most models for the origin of the primordial density per- 
turbation are based on a slowly-rolling, homogeneous infla- 
ton field. However, in view of our incomplete understanding 
of physics at the energy scales of inflation, combined with the 
increasing precision of cosmic microwave background (CMB) 
observations, it is important to consider other possibilities for 
the origin of the primordial density perturbation. Various al- 
ternatives to the orthodox scenario have been considered, in- 
cluding features in the inflaton potential |[ll|^[l,|4|,i5e 6], tran- 
sient effects of a pre-inflationary era 101, and trans-Planckian 
signatures |8]. In particular, a transient departure from slow- 
roll inflation has been considered, modelled by step feature 
in the inflaton potential lU IH [H 3] • Here we will consider 
a more general possibility: a transient departure from homo- 
geneity and isotropy in addition to slow-roll, modelled by an 
inflaton potential which is space-dependent during the step. 

The step feature considered in d 12. [H E'l was originally 
motivated by the idea that a second scalar field x|/ could cou- 
ple to the inflaton (|) and undergo a phase transition. The phase 
transition is triggered by the inflaton field in a manner sim- 
ilar to the transition of the waterfall field of hybrid inflation 
models ||^[l^. The i|/ field acquires a tachyonic mass and un- 
dergoes a phase transition towards the (|)-dependent minimum 
of its potential. The energy density as a function of (|) will then 
decrease until \(/ reaches its (|)-independent minimum. The ef- 
fect of this step down in energy density as a function of (|) is 
modelled by a step feature in the single-field effective poten- 
tial for (|). However, this reasoning also suggests that the step 
feature can have a space-dependence. Hybrid inflation-type 
pha se transitions commonly proceed via tachyonic preheating 
lUin . where tachyonic amplification of sub-horizon quantum 



fluctuations of \(/ results in an inhomogeneous energy density 
llTll[T2ll . In this case we expect the energy density as a func- 
tion of (|) to be space-dependent during the step. Therefore 
if a step feature can qualitatively model the effect of a phase 
transition on the inflaton, a space-dependent step feature is 
equally well-motivated. 

Regardless of any specific motivation, a space-dependent 
step feature can also be considered as a purely phenomeno- 
logical model for a transient breakdown of homogeneity and 
isotropy. The issue then is how the space-dependence generi- 
cally modifies the curvature perturbation relative to the homo- 
geneous step. 

Our paper is organised as follows. In Section II we intro- 
duce and discuss the space-dependent step model. In Sec- 
tion III we derive the equations for the classical and quantum 
modes. In Section IV we solve these equations for their effect 
on the curvature perturbation. In Section V we discuss our 
conclusions. 



II. A SPACE-DEPENDENT STEP FEATURE 

In order to study how space-dependence alters a homoge- 
neous step model, we will focus on the well-known (|)^ chaotic 
inflation step of y], [si |4|] . The potential of this model is given 
by 



l+c tanh 



(1) 



where b is the value of (|) at which the step occurs, d parame- 
terizes the duration of the step and c parameterizes the height 
of the step. We generalise this by adding a space-dependent 
term which modulates the step, 
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sin{kLx) sm{kLy) sm{kLz) , (2) 



2 



where 



1 



(3) 



The factor g{(^) parameterizes the magnitude and rate of dis- 
sipation of the inhomogeneities in the step, where (|)(f) is the 
classical background field. This corresponds to a lattice of in- 
homogeneities of characteristic length scale Xl = 2K/ki. Al- 
though approximately isotropic, a degree of anisotropy is in- 
evitable once inhomogeneity is introduced. In practice, in or- 
der to simplify the numerical analysis, we will consider an 
inhomogeneity which is planar in form, varying only in the 
ki direction. 



V( 
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-F((|))sin(kL.x) 



(4) 



In Appendix A we show that the mode equations for the 
planar inhomogeneity are similar to those of the approxi- 
mately isotropic inhomogeneity of Eq. (|2|. The second term 
in Eq. (HJl effectively replaces c with a space-dependent height 
for the step, c(l + g((|)) sin(ki.x)). We will consider the 
case where the inhomogeneities exist over the duration of the 
step and the variation of the height of the step with x is ap- 
proximately Cf, which can be modelled by setting g{(^) = 

(c/r/c)exp ^— ((|) — Z?)^ /c/^^ , where we assume Cf < c. This 

rapidly tends to zero for — > d. Note that we only use 
the background field in g{(^) since here (|) serves only to pa- 
rameterize the time-dependence of the inhomogeneities and is 
therefore not a quantum field. 



III. MODE EQUATIONS FOR THE CLASSICAL AND 
QUANTUM MODES 

There will be two contributions to the curvature perturba- 
tion from the space-dependent step. First, since the classical 
inflaton field is rolling in a space-dependent potential, it will 
acquire a classical spatial mode of initial wavenumber equal 
to the inhomogeneity wavenumber ki. Second, there will be 
quantum fluctuations of the inflaton amplified by the step as 
in the homogeneous case, but with additional features due to 
the space-dependence of the step. 

To study the combined metric and field perturbations, 
we introduce the Mukhanov variable u = —z1{_, where z — 
aisf/H and is the curvature perturbation fl3|. The space- 
dependence of the inflaton is defined to be entirely in the per- 
turbation m(x,t) about the homogeneous background. The ac- 
tion for u in the absence of the space-dependent potential is 
then UJ] 



^^2 



(5) 



where indices are with respect to conformal time T and co- 
moving coordinates x, and 



2H (j) 
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In this primes denote derivatives with respect to T and (|) = 
(|)(f). Note that (|) is treated as constant in the derivatives 
of Vstep but is subsequently set equal to (|)(f). Eq. ^ is of 
the same form as for a real scalar field in flat space with a 
time-dependent mass squared ml = — z /z- When the space- 
dependence is included, the effect is to modify the last term 



- — > — If' ^TTT sin (kL .x) — 2a^ sin (kL .x) m . (7) 
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We then derive the equations for the modes including space- 
dependence. Expanding m(x,t) in classical modes 



m(x,x) 



1 



(27r) 



3/2 



d^ku]i(T:)e 



ik.x 



(8) 



and substituting into the equation for m(x,t) with the space- 
dependence included, 



yd^F . „ , ^dF 



u — V^M= — u — a^ — ^sin(kL.x)M — — sin(kL.x) , (9) 



we obtain a system of coupled equations for the classical 
mode functions 

" 2 ^" '^^ '^^P I \ 

"k + K Mk Mk + TT —TTT ("k-ki — Mk+ki) 

z 2i d<Sf^ 



„3 

+ ^ ^ (27t)3/2 (5 (k - ki) - 5 (k + ki) ) = . (10) 
2i d^ 

The last term, which applies only to the modes with k — ik/,, 
is due to the rolling of the inflaton in the space-dependent po- 
tential, which creates a classical inflaton mode with k = ik/,. 

To quantize the modes, we will consider the action S with- 
out the linear term in u. The linear term only affects modes 
which are integer multiples of k/,. Therefore we can consider 
the modes with k ^ nk^ to be entirely quantum in origin and 
due to the terms in the action which are quadratic in m, while 
the modes with k = nk^ can be either quantum or classical, 
depending on which contribution is dominant. 

At times well before and after the step when the space de- 
pendence is effectively zero, the field may be expanded in 
terms of creation and annihilation operators with respect to 
the vacuum at an initial time, defined as A^; e-foldings before 
the end of inflation. 



m(x,t) 



d^k 



(27t 



3/2 



Mk(T)e'''-''flk + MkW- 



(11) 



The field operator m(x, t) and its canonical conjugate 7t(x, t) = 
M (x,t) satisfy the usual equal-time commutation relations and 
Hamiltonian field equations. 

The quantum modes for k 7^ nk/, satisfy Eq. ( fTOl i without 
the last term. 

" 2 Z a d F 

Mk+k Mk Mk+7r"JIT(«k-ki-Mk+ki) =0 . (12) 

Z 2.1 a'^^ 
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It is not immediately obvious that the quantum mode functions 
will satisfy the classical field equations. In the case without 
space-dependence, when substituting Eq. (fTTI) into the classi- 
cal field equation, the creation and annihilation operators fac- 
tor out. Including the space dependence introduces mode mix- 
ing and means this is no longer true. However, by perform- 
ing a unitary transformation to new canonical coordinates and 
conjugate momenta, it is shown in Appendix B that the late- 
time solutions for the quantum modes Mt are indeed the solu- 
tions of the classical mode equations Eq. ( fTSl l obtained using 
the quantum mode initial conditions. The initial conditions 
are obtained assuming a massless field with the modes of in- 
terest initially well within the horizon (k^ ^ a^H^), giving 



,-;(|k|x+9k) 



Mk -■ 



(13) 



where 9k is a random phase. There will also be quantum con- 
tributions to the modes at k = nk^. In this case we assume 
that the larger of the classical and quantum contributions dom- 
inates. 

We have included a phase 9k for each mode. In the case 
without space-dependence this phase has no physical effect 
and is usually set to zero. However, as we have shown, space- 
dependence implies that different modes couple, in which 
case the phases will affect the solutions for MkCt)- In ab- 
sence of phases, the quantum mode solutions Mk would gener- 
ate an ensemble of classical universes with the classical field 
mode probabilities determined from a Gaussian distribution 
with mean value |Mk| UM- We would then interpret the ob- 
served Universe as one of the ensemble of possible classical 
universes. When we include the initial phases and mode mix- 
ing at the step, each set / of initial random phases will generate 
a distinct solution for the modes at late times, Mk, /• Each of 
these solutions should be thought of as producing a distinct 
set of classical states which contributes to the total ensemble 
of classical states. Measurement of the state of the universe at 
late times will then select out a classical universe related to the 
solution Mk, , for one particular set of initial phases. Therefore 
in the following we will compute results for one set of random 
initial phases 9k. 



IV. NUMERICAL SOLUTION OF THE MODE 
EQUATIONS 

(i) Classical Modes 

The space-dependence of the step potential will imprint a 
classical space-dependence on the inflaton of wavenumber k^. 
The coupling of the modes in Eq. (fTOl l will then induce further 
classical modes of wavenumber nk^. In order to analyse these 
modes, we will expand the classical field as 



inki .X 



(14) 



which gives, on substituting into Eq. 

fl2 ^2p 
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Note that u'jj{i) shares the same dimensions with m(x,t) but 
not with the quantum modes Mk(T). We will consider through- 
out the case of a mode k = (A:, 0,0) in the x-direction, with 
ki — (A;/,, 0,0), where k and ki are positive. 

We have solved these equations numerically for the case of 
c ~ cp — 0.002, such that the spatial fluctuation of the step 
potential is of the order of the mean potential. (The maxi- 
mum value of c consistent with observations in the case of the 
homogeneous step is 0.003 jst].) 




FIG. 1: Im^'/^P as afunction of = for c = = 0.002 and /:t = 1, where 
k is in units of aH\sjEP- This coiresponds to a step with inhomogeneities of 
wavelength equal to the horizon at the step. 

In Figure [T] we show \u'jl/z\^ for the case of a step with 
inhomogeneity wavelength equal to the horizon at the step, 
ki = aH\^j^p. (We plot k in units of aH\^j^p.) A se- 
ries of classical modes are generated with wavenumber nki. 
However, the secondary modes at n = 2, 3, ... are very much 
smaller than the n — \ mode, which will dominate the cur- 
vature perturbation. The value of \u'jl /z\ for this mode is 
slightly larger that the conservative limit following from the 
observed value of < >, |m«V-z| < 10^"* (the relationship 
between the modes and the curvature perturbation is discussed 
in Appendix C), and therefore makes a larger contribution to 
< 'kJ' > than the homogeneous step. In this case the features 
in the curvature perturbation power spectrum will be quite 
different from the homogeneous step, with domination by a 
single perturbation mode of wavenumber k^ rather than the 
oscillations expected for the homogeneous step. 

For smaller values of cp we find that the height of |m','/zP 
decreases as c^. Therefore if the spatial fluctuation of the 
step potential is significantly smaller than the homogeneous 
step in the potential, the effect of the classical modes will be 
negligible. However, for Cf not too much smaller than c, the 



4 



curvature perturbation will be a combination of an oscillat- 
ing power spectrum of quantum origin plus a single classical 
mode at k = ki. 




FIG. 2: [w^'/zP as a function of A = nki for c = Cf = 0.002 and /t/, = 10. 
This corresponds to a step witli inhomogeneities of wavelength 10 times less 
than the horizon at the step. 
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|Mf,'/z| relative to the horizon-sized step. 

We conclude that the contribution of a classical spatial 
mode to the curvature perturbation is strongly dependent on 
its size relative to the horizon at the step. A significant contri- 
bution to the curvature perturbation is possible only if the size 
of the inhomogeneities is close to or larger than the horizon at 
the step. In this case we can expect to see a significant devi- 
ation from the smoothly oscillating power spectrum predicted 
by the homogeneous step feature. In particular, we expect the 
power spectrum will have a strong contribution from modes 
with wavenumber characterized by the inhomogeneity in the 
step. For the case of a less idealized step feature which is 
approximately isotropic and can be expanded as a range of 
modes with |k| ^ k^, the curvature perturbation will be dom- 
inated by a range of modes with |k| w rather than a single 
mode. 



(ii) Quantum Modes 

We next consider the purely quantum modes due to the 
quadratic terms in the action for u. It will be convenient to 
change variable to a = log(fl), in which case Eq. (fTZt be- 
comes 
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FIG. 3: |h;;'AP as a function ofk = nkt for c = cf = 0.002 and kt = 0.1. 
This corresponds to a step with inhomogeneities of wavelength 10 times 
greater than the horizon at the step. 

In Figures|2]and|3]we show the effect of varying the length 
scale of the inhomogeneities relative to the horizon at the step. 
Figure |2] shows |m','/zP for an inhomogeneity of wavelength 
equal to 0.1 times the horizon at the step. There is a dra- 
matic reduction in |m^//z| by a factor of 10^^ compared with 
the horizon-sized step. Since in the case of steps motivated 
by phase transitions we expect the inhomogeneities to be due 
to tachyonic amplification of sub-horizon quantum fluctua- 
tions, this suggests that for such steps the effect of space- 
dependence on the power spectrum will be negligible. This is 
an important conclusion, as it means that the oscillation fea- 
tures observed in jst] are a robust prediction of phase transi- 
tions during inflation, even if they are of the more complex 
tachyonic preheating type. In Figure [3] we consider the possi- 
bility that the space-dependence during the step is larger than 
the horizon. In this case we find a strong amplification of 



Ha = -47tG^^ 



(18) 



The subscript a denotes differentiation with respect to a. We 
have solved the system of equations Eq. ( fT6b , Eq. ( [TtI i and 
Eq. ( fTSl ) for a range of values of k of interest. To do this we 
have used a system of seven equations for each k, correspond- 
ing to full equations for Mk, Mk±ki and Mk±2ki, plus the equa- 
tions for Mk±3kL with the space dependence set to zero. Ran- 
dom values of the phase 9k were generated in order to study 
the effect of initial mode phases on the power spectrum. 

The initial conditions for the real and imaginary parts of Mk 
at k^ 3> a^H^ follow from Eq. ( fT3] l and are 



cos (9k 



IkiT) 



«k 



sin(9k+|k|T) 



-"ka 
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sin(9k 



k| 1 



2 e«// 



cos(9k 



|k|x) , 
(19) 

|k|T) , 
(20) 
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where Mk = + '"k- 

We solved for Mk, and used this to calculate |Mk/zp and the 
primordial curvature perturbation power spectrum, 



//2(k) = 




(21) 



Although the power spectrum in the case of an anisotropic per- 
turbation does not have a direct interpretation in terms of the 
mean squared curvature perturbation per logarithmic interval 
in k, we have plotted {k) with k = k^ in order to compare 
with the results for the homogeneous step [3]. As discussed in 
Appendix C, the more useful quantity for calculating the ob- 
served curvature perturbation is |Mk/zp, which we also plot. 
After the step has completed, the superhorizon modes grow 
proportional to z, so |Mk/zp is constant. It is in this limit that 
we calculate {k) and |Mk/zP for the superhorizon modes. 

The two components of Mk oscillate prior to leaving the 
horizon, with a frequency that increases as a^' as we go to 
earlier times. The background field (|) is initially slow-rolling 
and each mode is evolved from |k| >> aH to |k| << aH, 
The initial time for each mode is chosen differently in order 
to avoid the prohibitively small step size required for high 
frequency oscillations at very early times. A second-order 
Runge-Kutta algorithm was found to provide accurate solu- 
tions of the system of equations. 




FIG. 4: The curvature perturbation power spectrum for c = 0.002 and Cf = 0, 
where k is in units of aH\sTEP- Scales accessible to observations are between 
the dashed lines. 

Figure |4] shows the results for the case without space- 
dependence, Cf = 0, which reproduces the results of Jll]. 
In this and other figures we use c = 0.002, d — 0.01 and 
m = 2.5 X 10^^ (in units with li — c — G — I), where the 
value of m is determined by the COBE normalized value of 
. We consider the observable range of k to correspond to 
modes which leave the horizon from N = 60to N = 50, indi- 
cated by the dashed lines, with the step centered at = 55. 
The wavenumber which exits the horizon at the step occurs at 
the mid-point, A: = 1 in units of aH\sTEP- 

Figure |5] shows the curvature perturbation power spectrum 
rpj^ (k) with k — A:,, for the case of a space-dependent step 



PR(k) 10-8 




FIG. 5 : The curvature perturbation power spectrum for inhomogeneities of 
wavelength equal to the horizon at the step, kt = \, for the case Cf = c = 
0.002. Resonances are observed at k = k/. and k = 2ki. The apparent broad- 
ening of the power spectrum is due to the random initial phase for each mode 
combined with mode coupling. 




FIG. 6: |«k/zp for inhomogeneities of wavelength equal to the horizon at 
the step, Ul = 1, for the case Cf = c = 0.002. 



of wavelength equal to the horizon at the step (Al = 1). In 
Figure |6] we show the corresponding plot of |Mk/zp. There 
are two notable new features in the power spectrum. In ad- 
dition to the oscillation feature of Figure IH we observe (a) 
sharp resonance-like features and (b) an apparent broadening 
of the spectrum relative to the line observed in the homoge- 
neous case. 

(a) Resonances: The resonance-like features are found to oc- 
cur at k = ki, k = 2ki and k = 3kL- We expect that had we 
used modes up to Mk+nk^ in the calculation of Mk, we would 
observe n resonances, decreasing in magnitude as n increases. 
While having a large amplitude, the width of the resonances 
in k is very narrow, which greatly diminishes their effect on 
the curvature perturbation. 

The origin of the resonance is the coupling of a mode Mk to a 
much larger superhorizon mode Mk-k/,, where |k — ki| ^ aH. 
Since the superhorizon growing mode has been strongly am- 
plified, with Mk-„ki °^ a since exiting the horizon, it pro- 
vides a strong driving term in the Mk mode equation, produc- 
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FIG. 7: A close-up of the k = resonance of Figure|6] 



ing rapid growth of the mode. The longer the wavelength 
of the k — ki mode, the larger will be the amplification of 
Mk-ki- Therefore the Mk resonance will tend to large values 
as |k — ki| ^ 0. A close-up of the k = resonance of Fig- 
ure |5]is shown in Figure |7] (A lower bound 5 x 10^^ has 
been imposed on |k — k^j when plotting the figures.) We find 
from Figure |7] that |Mk/zP °= |k — kij^-'. The resonance de- 
pends only on |k — k^j because the amplified superhorizon 
mode driving the growth of the resonance depends only on 
|k-ki|. 

Although resonance due to coupling to large superhorizon 
modes is a theoretically interesting consequence of the space- 
dependent step, we believe it will make a negligible contri- 
bution to the curvature perturbation compared to the classical 
modes with k = nk/,. In Appendix C it is shown that a quan- 
tum resonance produces an effective classical mode given by 



the horizon at the onset of inflation. In this case the logarithm 
in Eq. ( |23] ) can only enhance the prefactor by a few orders of 
magnitude at most. If s were significantly larger than 3 then a 
large enhancement might be possible as bk 0. However, we 
find no evidence for s > 3 from our numerical results. 

From Figure |2l assuming that \u]i/z\^ °^ k ^ for all k > bk, 
we find A3 — 4 X 10" '"^ and so for « = 1 



:2x lO^^^ln 



5k 



(24) 



For any reasonable value of 5k this will be much smaller than 
the « = 1 classical mode and also much smaller than the ob- 
servational upper bound of 10"^. 




FIG. 8: |«k/z|" averaged over many sets of random phases for inhomo- 
geneities of wavelength equal to the horizon at the step, kt = i, for the case 
Cf = c = 0.002. The apparent broadening in Figure[6]is no longer seen. 
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|k-nkL| . (22) 



The dependence of |Mk/zp on Ic is critical. If we integrate 
Eq. (l22l) from a lower bound k = 8k to ki ^ bk, then with 
|Mk/zP —As/k^ we obtain 
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Therefore if 5 = 3, as our numerical results indicate, the log- 
arithmic dependence on the cut-off will prevent a large en- 
hancement as 5^ 0. We generally expect a physical lower 
bound on dk to exist, corresponding to the largest unsup- 
pressed superhorizon mode. For example, the amplitude of 
superhorizon modes in (|)^ chaotic inflation will become sup- 
pressed by effect of the finite mass term, while more generally 
the largest amplified mode will correspond to the mode exiting 



FIG. 9: Close-up of Figure[6] showing the discrete random fluctuation of the 
power spectrum with each step in k. 

(b) Broadening: The apparent broadening of the power spec- 
trum is due to the random initial phases of the modes. For each 
value of k, the random phase produces a value of iP,^ '''^(k) 
which is shifted either above or below the homogeneous re- 
sult. In Figure [8] we show the same plot as in Figure |5] 
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but averaged over many sets of initial random phases. Seen 
close up, the power spectrum is actually a series of discrete 
points, randomly above or below the mean spectrum, at val- 
ues of k separated by the step size of the numerical solu- 
tion. This can be seen in Figure |9] (a close up of Figure |6] 
for 1 .75 X 10"^ < < 2.05 X 10"^). There is no apparent cor- 
relation between for different values of ^ - it is a random 
scatter about the result for zero phase. This scatter is due to 
mode coupling, with each mode having a completely random 
initial phase. In the case of a homogeneous step there is no 
mode coupling, therefore the initial phases can be factored 
out of the mode equations and have no observable effect. We 
have produced solutions using several different sets of random 
initial phases and find no qualitative difference in the result- 
ing power spectrum, although the detailed pattern of random 
fluctuations depends on the initial phases. 

In practice, when translating the primordial power spectrum 
into physical observables, we average over the contributions 
of modes for a finite range of k. Since the random fluctuations 
will average out to zero over this finite range of k, it is likely 
that the observable effects will be indistinguishable from the 
case where all initial phases are set to zero (which is equiv- 
alent to averaging over many sets of initial random phases). 
However, the transient coupling of modes during the step may 
also lead to some correlation of the phases of the modes. This 
could generate non-Gaussanity in the curvature perturbation. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have been concerned with how space- 
dependence of the inflaton potential can manifest itself in the 
curvature perturbation. 

A space-dependent step feature provides a model for a tran- 
sient breakdown of slow-roll, homogeneity and isotropy dur- 
ing inflation. We have studied how the curvature perturbation 
of the well-known (|)^ chaotic inflation step feature of JK] is 
modified by space-dependence. 

We find that the effect of space-dependence depends 
strongly on the size of the inhomogeneity in the step potential 
relative to the horizon at the step. The curvature perturbation 
can be dominated by the effect of the space-dependence if the 
inhomogeneities are close to or larger than the horizon, but 
their effect is strongly suppressed if the inhomogeneities are 
small compared with the horizon. 

One motivation for step features is a phase transition oc- 
curring during inflation. If the phase transition was of the 
tachyonic preheating type, we would expect the step in the in- 
flaton potential to be strongly space-dependent. However, the 
space-dependence would typically be on sub-horizon scales, 
coming from tachyonic amplification of sub-horizon quantum 
fluctuations. If the space-dependence is much smaller than 
the horizon, our results imply that space-dependence will have 
a negligible effect on the curvature perturbation power spec- 
trum. Therefore the smoothly oscillating power spectrum ex- 
pected for the homogeneous step is likely to be a robust pre- 
diction, typically unaffected by space-dependence during the 
phase transition. 



We also considered the possibility that the space- 
dependence is on scales equal to or greater than the hori- 
zon at the step. In this case the curvature perturbation can 
be dominated by the effect of space-dependence. For the 
simple model of an anisotropic step with inhomogeneities of 
wavenumber ki, the curvature perturbation will have a sin- 
gle classical mode of wavenumber ki superimposed on the 
oscillating power spectrum from the amplification of quan- 
tum fluctuations during the step. As a result, the signature 
of space-dependent step features can be quite different from 
the smoothly oscillating power spectrum expected for a ho- 
mogeneous step, with the curvature perturbation dominated 
by modes characterized by the wavenumber of the inhomo- 
geneity. 



The quantum modes are also modified by the effect of 
space-dependence and mode coupling. In particular, the quan- 
tum curvature perturbation has resonance-like features at in- 
teger multiples of the inhomogeneity wavenumber ki. These 
are due to coupling to strongly amplified superhorizon modes. 
Although the resonances appear large in a plot of the curva- 
ture perturbation power spectrum, they make only a very small 
contribution to the curvature perturbation compared with that 
from the classical modes, as a result of their narrow width. A 
second effect of mode coupling during the step is that the cur- 
vature perturbation becomes dependent upon the initial phases 
of the quantum modes, producing a scattering ('broadening') 
of the curvature perturbation modes as a function of k. How- 
ever, this random scattering is likely to average to zero in any 
physical observable. 



There may be other effects of a space-dependent step. One 
is that the transient coupling of modes during the step could 
induce some correlation between the phases of the modes, 
generating a degree of non-Gaussianity in the curvature per- 
turbation. It has been shown that non-Gaussianity is gener- 
ally produced by step features in single-field inflation models 
lllSlllql even in the absence of space-dependence. In addition, 
an oscillating potential may produce non-Gaussianity in sub- 
horizon modes at the step 1 16]. It will therefore be important 
to establish if non-Gaussanity from space-dependence can be 
distinguished from these other sources. A second effect is that 
an anisotropic space-dependent step feature could provide a 
mechanism to generate primordial anisotropy in the curvature 
perturbation. It has recently been claimed that there is a large 
scale statistical anisotropy in the observed CMB [ITtIi. One 
explanation of this is the existence of a large superhorizon 
mode in the curvature perturbation along a particular direc- 
tion k ifTsIl . This would require a mechanism to break isotropy 
and to amplify a single mode. Generation of an anisotropic 
space-dependent step of wavenumber k^ would be one way to 
introduce such a feature into the curvature perturbation, with 
modes with k w k/, being generated by the step. Although 
an inflaton field cannot account for the observed anisotropy 
without introducing too much CMB inhomogeneity, a slowly 
rolling curvaton field with a space-dependent step in its poten- 
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tial might be able to ifTsI] ' . A method of checking the pres- 
ence of large scale anomalies in the observed CMB tempera- 
ture has been suggested in lEoll . The polarization-temperature 
correlations presented there could be a way to test the step 
feature presented in this paper. 
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In this Appendix we prove that the classical field equation 
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Appendix A: Mode equations for an approximately isotropic 
distribution of inhomogeneities 

In our analysis we have considered mode equations for an 
anisotropic planar inhomogeneity in the step potential. In this 
appendix we generalise the mode equations to the case of an 
approximately isotropic distribution of inhomogeneities. In 
this case the step potential is given by 



V((\>)=V,s 



step V 



F{(^)sm{kLx)sm{kLy)sm{kLz) . (A-1) 



This corresponds to a 3-D 'lattice' of inhomogeneities of di- 
ameter 2k /Xl. Strictly speaking this also violates isotropy, 
but far more weakly than the planar inhomogeneity. Follow- 
ing the same steps as before, the classical mode equations are 
now given by 



1 2 Z 

k Mk Uk 



d^F r 



(2/)3 t/(^2 



k+kO 



"k-kV ~ "k+k! 



dF 



'k-k 



^k+k 



"k-kj^ 



«k+kv 



{2iy £/(j) 



(§k,kO 



--k.-k' 



'k,k; 



(5k,ki 



5k,-k! 



Jk,-k-; 



, (A-2) 



where 



= {kL,kL,kL), 



{—kL^ki^—ki) and k|^ = {—kL,—kL,ki). The qualitative 
behaviour of the solutions of Eq. ( IA-2b is expected to be 
similar to the case of the planar inhomogeneity described by 
Eq. ( fTOl ). with the difference that each mode now directly 
couples to eight other modes rather than two. In particular, 
we expect that resonant behaviour will still occur for any 
mode k such that |k± k^| is small for / = 0, x,y or z. 



' The possibility of a curvaton with transient breakdown of slow-roll has 
recently been explored in Il9il . 



1 



d'^x 



(V«) 



2 2 

- m„u - 



■m\ sin(kL.x)M^ 



-.dF 

-2a -rrsin(kL.x)M 



diSf 



(B-1) 



where = —z jz and = cp-d^F j di^ . The term linear 
in M drives the growth of the modes with k — ik^ and the 
modes with k = nk^ which are coupled to these. The growth 
of these modes, which is due to the classical evolution of the 
inflaton field in the space-dependent potential of the step, may 
be treated classically. We will therefore quantize the action 
without the linear term, and consider the classical evolution 
of the modes with k = nk/, separately. It will be convenient to 
use discrete modes in a box of volume V, 



m(x,t) 



7i:(x,t) 



1 



(k.x 



(B-2) 



(B-3) 



where 7i = m is the conjugate momentum of u and k, t) = 
<7^(k, t). The Hamiltonian for the modes in the absence of the 
linear term is then 

^ = f/''(k,T)/^(k,T) + (k2 +m2)^+(k,T)^(k,T) 



-^^^(k-Kk^,T)^(k,T) - ^^^^(k,T)^(k- 
LI LI 



This can be written in the form 
1 



(B-4) 



(B-5) 



where P and Q are column matrices with entries /?(k,T) and 
^(k, t) labeled by k and M is a Hermitian matrix. The terms in 
Eq. ( IB-41 ) proportional to m\ are due to the space dependence. 
Since these terms mix modes with different k, we cannot di- 
rectly quantize the Hamiltonian by expanding m(x,t) in terms 
of creation and annihilation operators with mode functions Mk 
which satisfy the classical field equations. Instead, we first 
transform p and q to new canonical coordinates and momenta 
which allow the Hamiltonian to be expressed as a sum over 
uncoupled quadratic Hamiltonians for each k. To do this we 
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make a unitary transformation of q and p to new variables s 
and ps. 



5=f/e^.9(k,T) = c/. .,^(k,T) 



P,^[/*P^;,,(k,T)=t/*^,/,(k,T) 



(B-6) 



(B-7) 



where [/ is a unitary matrix with rows and columns labeled by 
k which diagonalizes M. The Hamiltonian can then be written 
in the form 



1 



H^^Y. pi(k,T)/,,(k,T) + a)^/(k,T)^(k,x) 



(B-8) 



where a)^(k,k£,T) are the eigenvalues of the matrix M. In 
addition, s and ps are conjugate variables, as can be seen from 
the conjugate nature of q and p, 

\ps (k, T) , ^ (1, T)] = f/„, b(k' , T) , ^(l' , T)] 



[p,(k,T),i(l,T)] =-/5i.k 



Since Eq. ( |B-8l l is quadratic in the modes, we can quan- 
tize the system by expressing p^ and s in terms of annihila- 
tion operators and mode functions ^(k,!) = Ui^{i)d^, where 
Mk(x) will satisfy the classical equation for i(k,T), obtained 
using the Hamiltonian form of the equations of motion for 
canonical coordinates and conjugate momenta, p^ — —dH/ds, 
s = dH/dp,: 



/(k,T) = -a)^i(k,T) 



(B-9) 



At times well before and after the step, U]i{x) and Mk('c) are 
the same, since then s(k,T) ^ q{k,z). Therefore the quantum 
initial conditions for Mk(T) are the same as those for MkCt), 
given by Eq. ( fT3] l. By rotating back from S to Q via an inverse 
unitary transformation, the classical equations Eq. ( |B-9t be- 
come the classical equations Eq. ( fTOl i. Therefore the solution 
of Eq. (|B-9t for quantum initial conditions is the same as the 
solution of Eq. (fTOl i with quantum initial conditions given by 
Eq. (fT3] l. At times well after the step has occurred, Mklt) may 
be interpreted as the quantum mode function with the field 
given by Eq. (fTTT l. since theni(k,T) q{k,z). 



Appendix C: Relating Modes to tlie Curvature Perturbation 

(i) Power Spectrum for Anisotropic Perturbations 

The power spectrum of m(x,t) may be defined in terms of 
the spatial average of m(x,x)^. Setting the initial phases to 
zero for the moment, from Eq. ( fTTI ) this is given by 



(m(x,t)^) = {0\u{x,zy\0) = 



(2^ 



\Uk\ 



(C-1) 



In this we are applying the ergodic principle that classical spa- 
tial averages over a member of an ensemble are equal to quan- 
tum ensemble averages at fixed x. The mean value of m(x,t)^ 
is related to the power spectrum by 



(«(X,T)2) 



g»(k) 
471: IkP 



(C-2) 



Since in most cases the power spectrum is independent of the 
direction of k, iPuik) is usually written as a function of A; = |k|. 
In the anisotropic case we can use the more general definition 
of Eq. ( IC-21 ). Comparing with Eq. iC-li we then obtain 



rP^Ck) 



27l2 



«k 



Therefore, since ^ ~ —u/z, we find 



^ g»(k) _ |kp 

z2 27t2 



(C-3) 



(C-4) 



In our calculation of the power spectrum, we apply 
Eq. dC-ll ) to the solution Mt corresponding to a particular set 
of random initial phases. The average of Eq. (IC-ll l is then an 
ensemble average of u{x,z)^ at fixed x over the classical uni- 
verses generated by u^. Since the observed power spectrum 
is assumed to correspond to that obtained in one member of 
the ensemble by spatially averaging m(x,t)2, we are in prac- 
tice applying the ergodic principle to the ensemble generated 
by each set of random phases i.e. we first fix the initial phases 
and then apply the ergodic principle to the resulting Mk('c)- 



(ii) Relating classical and quantum modes to the curvature 
perturbation 

In the case of isotropic fluctuations, the power spectrum can 
directly interpreted as the contribution to < > per loga- 
rithmic interval in k. However, we are considering here an 
anisotropic step characterised by the direction k^. In order to 
interpret in this case, we can discretize m(x,t) into classical 
modes 



k, 



In this case 



c/|2 



For quantum modes we have 

(«(x,T)2)^(0|«(x,T)2|0)=|-^|«k|' 



(C-5) 



(C-6) 



(C-7) 



To relate the quantum resonant modes to classical modes we 
can write this as a sum of the form Eq. iC-6i with 



,.cl I' 
Steffi 



1 



(27t)^ J-kJl 



kL/2 



ukl^d^ki ; k,=k-k;. (C-8) 
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This breaks the integral Eq. iC-7\ into cubes of side cen- 
tered on points k, = {a,b, c)ki with a, b, c integers. In so far as 
the integral Eq. ( IC-7b is dominated by the contributions of the 
resonances, we can then identify \u^, | in Eq. ( IC-81 ) with the 
amplitude of effective classical modes in the sum Eq. dC-Sb 
with k; = nkt. 

The resonances are a function of |k — «ki|. Therefore the 
integral over the volume of the cube Eq. ( IC-81 ) can be ap- 
proximated by an integral over a sphere of radius k^. With 



,,cl 



and k = k — nk^, Eq. dC-Sb becomes 



\ui.\^{k)4nk ^dk 



(C-9) 



Since we can choose k along any direction, we can assume 
k = kx — nki with kx > nki- This allows us to directly compare 
the quantum modes with the classical modes. 

The curvature perturbation due to the classical modes is 
then 



k„ 



(k„ .X 



(C-10) 



of wavenumber nk/,, which will manifest itself as anisotropic 
density perturbations in planes normal to the k^ direction. 

To constrain the amplitude of the classical modes m^', we 
can interpret the observational constraint on Pj^ as a constraint 
on < 2^^ >. For the case of isotropic perturbations. 



< !K. > - 



dk 



(C-12) 



where k,„i„ and ^,„(„ correspond to the largest and smallest 
scales observed in the CMB. Using the COBE normalized 

value, P^^{k) — 4.8 x 10^^ (essentially independent of k for 
a near scale-invariant spectrum) and \n{k,nax / kmin) ~ 10, we 
obtain 



Since the contribution from the classical modes is 



(C-13) 



(C-14) 



where 



(C-11) 



Therefore the planar inhomogeneity in the direction k^ will 
result in a curvature perturbation which is a sum of modes 



we will therefore impose the conservative constraint 



< 10-4 



(C-15) 
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